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Abstract 

We show that any sphericahy symmetric spacetime locally admits a maximal 
spacelike slicing and we give a procedure allowing its construction. The construc- 
tion procedure that we have designed is based on purely geometrical arguments and, 
in practice, leads to solve a decoupled system of first order quasi-linear partial differ- 
ential equations. We have explicitly built up maximal foliations in Minkowski and 
Priedmann spacetimes. Our approach admits further generalizations and efficient 
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1 Introduction 



In Lorentzian geometry, a maximal hypersurface is one that is spacelike and has vanishing 
mean extrinsic curvature, K = 0, K being the trace of the extrinsic curvature of the 
hypersurface. The name comes from the fact that the induced volume functional reaches 
a local maximum with respect the variations that keep fixed a given boundary. Maximal 
hypersurfaces were considered by Lichnerowicz ^ to solve Einstein's constraint equations, 
giving motivation for subsequent studies on the subject .^^^^^^^^^^^^ In fact, the existence 
of maximal hypersurfaces is extensively used in Mathematical Relativity.^ This property 
is a very simple geometric assumption allowing one to establish general results for a very 
large variety of spacetimes, for instance, local or asymptotically stationary or conformally 
flat spacetimes. 

In this paper, we will use the term maximal slicing when referring to a (non intersect- 
ing) family of spacelike maximal hypersurfaces which locally foliates a certain domain of 
spacetime. This type of slicing has very nice properties as, for example: i) the well-known 
singularity avoidance capabihty,!^ ii) it is well adapted to the propagation of gravitational 
waves, ^^^^^ and, iii) it gives the natural Newtonian analogous when, in addition, con- 
formal flatness is imposed on each slice Maximal slicing condition has been recently 
used in the Fully Constrained Formulation of Einstein equations derived by the Meudon 
group,^I2IIll in the field of Numerical RelativityJ ^^ ^ ^^ I The study of maximal hypersurfaces 
is not only interesting to carry out practical applications, but also, to solve fundamental 
questions in Relativity. Consider, for example, the proof of the positiveness of the mass 
for isolated gravitating systems, that was firstly proved by assuming the existence of a 
maximal spacelike hypersurface,^^ and only then, it was finally established by removing 
this reasonable hypothesis.^^ 

The study of maximal slicings has been mainly motivated by the resolution of the 
initial value problem.I^^ This idea was popularized by Smarr and York.^^HSI Then, for a 
given class of spacetimes, the construction of a family of slices with a given property 
(e.g., maximal slicing and/or conformal flatness) may be dynamically accomplished in 
an evolution scheme: assuming the existence of an initial 3-surface with that property, 
Einstein equations are used to evolve locally the corresponding additional constraint. 

This strategy can be applied both analytically and numerically, as it has been done by 
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Estabroook et al.^ and by Beig and MurchadhaPH to construct maximal slicings in the 

extended Schwarzschild geometry. This maximal foliation of the Kruskal-Szekeres black 

hole was obtained independently by ReinharlP^ from a different approach. Reinhart's 
99 

procedure does not make any use of Einstein equations, being essentially a geometrical 
(non dynamical) approach. 

Similar examples of maximal slicings have been examined by MalecP^ in the Reissner- 



Nordstrom geometry, and by Petrich et al^^ in dynamical spherically symmetric space- 
times (SSSTs) of dust collapse scenarios, by using isotropic coordinates. Eardley and 
Smari^^ obtained general results and presented some conjectures about the existence of 
Cauchy time functions whose level surfaces are maximal (or foliations with constant mean 
extrinsic curvature) with direct applications to the Lemaitre-Tolman-Bondi geometry. Of 
course, the study of the existence (and uniqueness) of other types of foliations deserves a 
lot of interest in SSSTs and more general spacetimes. That is the case of synchronizations 
by flat instants (Painleve-Gullstrand slicings) and their generalizations.l^ ^ l ^''^l^^^^^ l "^^ l "^^ l "^ ^ 
Recently, foliations with constant mean curvature have been constructed by MalecP^ in 
the Schwarzschild geometry. 

In spite of their extended use, the existence of maximal slicings in SSSTs has been only 
established for vacuum and for some particular energy contentsl ^^ l ^^ l ^^ l ^^ l ^"^ ^ ^^ ^^ There 
is, as far as we know, no theorem stating that always it is possible to build a maximal 
slicing in a SSST. 

In this work, we aim to prove the local existence of maximal slicings in an arbitrary 
SSST. We will follow a purely geometrical approach, independent of Einstein equations, 
complementary to the standard time evolution strategyj ^'^ ^ ^^ ^ '^^ ^ '^'^ ^ '^^ l '^^ l Although our 
study is independent of the field equations, one by-product of our approach, which could 
be of interest in the field of Numerical Relativity, is that it provides a means to assess 
complex and sophisticated 3D numerical codes built to solve Einstein equations. Indeed, 
most of the current 3D present-day Numerical Relativity codes are written in Cartesian 
coordinates. Two classical, standard test-beds in spherical symmetry are the long-term 
stable evolution of a Schwarzschild black hole (vacuum), and the Oppenheimer-Snyder 
coUapsd^of dust (matter). In both cases the builder of a Numerical Relativity code aims 
at the recovering of the corresponding well-known solutions in spherical symmetry. This is 
an obvious assessment (in the 3D case) of a code. Let us consider the Oppenheimer-Snyder 
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analysis of the collapse of a ball of dust; it is the simplest approach to the current studies 
on the complex phenomena of the stellar core collapse (in the astrophysical scenario of 
hydrodynamical supernovae). Hence, for a given 3D numerical code, written in Cartesian 
coordinates, the Oppenheimer-Snyder solution is a nice test-bed to asses the performances 
of that code, but also to calibrate it against any spurious gravitational radiation released 
as a consequence of grid effects. Since our approach provides a procedure to generate 
a maximal slicing in a SSST, this, in practice, means that any SSST, not necessarily 
satisfying the maximal slicing condition, can be reformulated in terms of this slicing and 
used to check a numerical 3D solution generated by a code written in Cartesian coordinates 
or not and satisfying the maximal slicing gauge. 

Our analysis has focused on SSSTs in standard coordinates. An intrinsic characteri- 
zation of SSSTs in terms of metric tensorial concomitants has been very recently accom- 
plished by Ferrando and Saez.l^Sl The study of SSSTs is of particular interest in different 
domains of general relativity: classification of exact solutions of Einstein equations, sim- 
plified models of astrophysical and cosmological scenarios,^^ test-bed solutions for fully 
general-relativistic time-dependent numerical codes which evolve matter in strong grav- 
itational fields.^^^l^^I^S12S] SSSTs provide the background for a perturbative theory. All 
these reasons advise us to start our analysis in such a simple symmetry and to defer in a 
future paper its extension to less strongly symmetric spacetimes. 

The paper is organized as follows. In Sec. [2] we consider (in a SSST) vorticity-free and 
non-expanding observers, that is, those that are orthogonal to maximal slicings. In Sec. [3] 
we analyze the existence of a coordinate transformation allowing a new coordinate time 
that has associated a maximal slicing. In Sec. H] we derive a decoupled system of three 
first order partial differential equations that proves the local existence of a maximal slicing 
and, in Sec. we provide a general procedure allowing its construction. This procedure 
is illustrated in Sec. [6] considering the most simple geometry: a flat spacetime, solving 
the problem analytically. The procedure is also illustrated in Sec. [7] considering a closed 
Friedmann spacetime, which is used in physical applications of dust collapse. Finally, in 
Sec. [HI we summarize the results. It will be understood that the title of the sections always 
concern spacetimes with spherical symmetry. A short communicatioiJ^ of preliminary 
results of this work was presented at the Spanish Relativity meeting ERE-2009. 
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2 Kinematical approach to maximal slicings 



Without loss of generality, we start with the canonical form of the metric of a SSST: 

ds^ = Adf + 2Cdtdr + Bdr^ + DdQ^, (1) 

(see, e.g. the book by Plebahski and KrasihskPHl) where dfl^ = d^^+sin^ Odip'^ is the metric 

of the 2-sphere, and A, B, C, D are smooth functions of t and r such that AB — < 0, 

which guarantees the Lorentzian character of the metric. In addition we choose the 

signature (— , +, +, +), and accordingly D > Q. Partial derivatives with respect to r and 
df df 

t will be denoted as — — = /' and — - = / respectively. The spatial metric 7i, induced on 
or ot 

the hypersurfaces Sj, defined hy t = constant, is •jij = diag(i?, D, D sin^ 6^), where 5 > 
since we are considering spacehke hypersurfaces. Let n = n^dfj, be the future pointing 
timelike unit normal to the hypersurfaces Sj, 



1 f d C d\ , 



The mean extrinsic curvature K of is 



DVC^ - AB ' 



(3) 



where V is the covariant derivative with respect to the spacetime metric given by Eq. ([T]), 
and /i G {t,r,6,ip}. Simplifying this relation, we obtain 

That the surfaces t = constant are maximal means that the maximal slicing condition 
i^' = is fulfilled. In the following, we assume that A,B,C and D are known functions. 
We aim to prove the local existence of a maximal slicing in spherical symmetry. From 
a kinematical point of view, we look for the existence of a non-expanding observer u on 
the 2-surfaces defined by constant 9 and ip, i.e., a unit and future directed radial timehke 
field u = = {u^{t, r),u^(t, r), 0, 0), such that V^u^ = or, equivalently. 



*DVC2 - AB) ' + {u'DVC^ - AB)' = 0. (5) 



The condition g{u,u) = —1, with respect to the spacetime metric g^u given by Eq. ([T]), 
leads to the following relation between the components of u: 

A(u*)2 + 5K)2 + 2CmV = -1. (6) 



Eq. (|6]) allows to express u"^ in terms of -u* and the metric components when the associated 
discriminant is not negative, that is, 

Substituting the expression for in Eq. (E]), one obtains a first order partial differential 
equation for m*. In order to guarantee the existence of a maximal slicing, the solutions of 
this equation have to satisfy the above algebraic constraint ([7]) on u*. 

In the next section, we introduce a change of coordinates |t = i{t, r), f = f{t, r), 9, ip} 
and we require that the new coordinate hypersurfaces t = constant be maximal. As we 
show below, no additional algebraic constraint must be considered when the conditions 
assuring the existence of this coordinate transformation are properly taken into account. 



3 Main equations of maximal slicings 

According to the previous comment, we need to introduce two fields, X and Y, such that 
the commutation relation [X, F] = be satisfied. This condition assures the existence of 
two coordinate parameters, namely t = t{t,r) and f = r{t,r), such that 

X = l,Y = l. (8) 
dt dr ^ ' 

We decompose these fields 

Y = XY, = 1, (9) 

X = aY + bY^, F-F^ = 0, (10) 

d 

with 6 7^ and A > 0, and fixing equal to one the coefficient of — in the decomposition 

dt 

d d 

of y-*- in terms of — and — . The condition = 1 is equivalent to 
ot or 



fill 



being / an unknown function to be determined, e = ±1 and = —AB + > 0. 
Moreover, Y ■ Y-^ = leads to 

r^ = |- + g^, Q = B-'(-c + ^ML=]. (12) 

dt ^dr' ^ \ ^pp j^B 
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Consequently, the resulting fields are 



X = iaf + b)^^ + iaP + bQ)^, 



(13) 
(14) 



where we have taken into account Eqs. ([2]) and ffTTl) . Using this decomposition, we will 
obtain an equation, involving only the unknown /, that gives a maximal family of hyper- 
surfaces t = constant. Notice that Eq. f|T^ makes clear the meaning of e: relative to the 
Eulerian observer n, the radial component of Y is outgoing or ingoing according to e = 1 
or e = —1. 

The Jacobian matrix of the above change of coordinates can be written, in terms of 
the previous decompositions, as 



i i' 



/ dt dt \ 

dt df 

dr dr 

\ dt df I 



(bpX) 



-1 



-1 

' af + b Xf 
aP + bQ XP 

XP -Xf 
-{aP + bQ) (af + b) 



(15) 



where p = P — fQ 



7^ 0. In order to ensure Eq. ([H]), we must require that 



[X, y] = 0, and this condition is expressed according to the following result. 

Lemma 1: For the fields X,Y given by Eqs. fT^) and [14\ ), the condition [X,Y] = 
is equivalent to 



(16) 
(17) 



"p" 




[-/] 


bp 




bp 



aP + bQ' 


1 


'af + b 


bp X 




bp X 



To avoid confusion, we denote with K the trace of the extrinsic curvature of the new 
hypersurfaces t = constant. From Eq. (jl]), the maximal slicing condition for the new 
hypersurfaces, K = 0, is equivalent to 



X(Y') + 2Y'^^ 



2Y{X ■ Y) + 



Y(Y^ 



D 



X ■ r = 0. 



'1^ 



Taking into account Eq. f fTSj) and the commutation relation, we have 3 equations for 4 
unknown functions, a, 6, /, A. But, due to D transforms as a scalar under the considered 



coordinate transformation, we can add without loss of generality, the coordinate condition 

f^Y^ = D, (19) 



saying that the metric on the hypersurfaces t = constant is written in isotropic confor- 
mally flat form. Notice that Eq. (1191) guarantees that the congruence of coordinate lines 

associated with the field X = — has minimal distortionj ^ ^ ^^ l 

dt 

From Eq. f|T8|) . taking into account the conformally flat coordinate condition (|T9|) . we 
obtain the following result. 

Lemma 2: The condition that the hypersurfaces t = constant are maximal (K = is 
equivalent to 



2Y {X-Y)- 3X (F^) + 



4 Y{Y 



2\ 



- + 



X • y = 0, (20) 



/ y2 

where X, Y and f are defined by Eqs. ( [7^ and ( flgj) . 

Lemma 2 expresses the maximal slicing condition in terms of the fields X, Y . These 
fields have to satisfy the commutation relation given by Lemma 1. In the next section, 
we rewrite both conditions in terms of the decompositions (fT3|l and (JH]). 



4 Local existence of maximal slicings 

In this section, we deduce the theorem stating the local existence of maximal slicings 
in SSST. The proof consists in reducing the required conditions into a set of first order 
partial differential equations. 

In terms of the decompositions (fT3l) and ( fT4l) . Eqs. ( IT9l) and (l20l) are expressed as 



A = ^D/f (21) 



Eq. (121]) can be viewed as a definition of f in terms of A, and from this equation we can 
derive 

Y ^ _ af + b A = A + <^P + ^Q ,23) 
A 2D bpy/D' A 2D bp^fD ' ^ ' 

where we have taken into account the expression of the Jacobian matrix given by Eq. fllSp . 

At this point we have 5 equations, i.e., Eqs. f lT6|) . f|T7|) . f l22|) and fl23|) . for 4 unknown 
functions, i.e., (a, 6, /, A). Eqs. (|23|) are equivalent to give quantity a in terms of the other 
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unknowns, 

^ / i -^ o\' h -i-nn'\ 

(24) 



^ ^ / A + QA' D + QD' 



and to the following equation 

Once a has been eliminated, it can be seen that Eq. f fT7|) can be obtained from Eq. fl25|) . 
From Eq. (^^, the field X given by Eq. f[T^ is expressed as 

X = .pyD(-A'| + A^), (26) 

where A = log A/V^. So we have reduced the problem to solve 3 equations, i.e., Eqs. f|T6|) . 
( l22l) (replacing a), and ( l25i) . for 3 unknown functions, i.e., (6,/, A). Eq. ( |T6i) can be 
rewritten as 

f\ + P^l = P'-pP- + f-fl. (27) 
p p 

As it will be shown below in Eq. fl25]) . / can be obtained from an equation which does not 
contain other variables. Then, Eq. (127|) is a first order partial differential equation (PDE) 
for b. 

Eliminating a and its derivatives with Eq. ( |24l) . and making use of Eq. ( l27l) to eliminate 
b, Eq. ( I22l) is reduced to a second order PDE involving A and /. Finally, making use of 
Eq. (j25|) for, firstly, reducing the order of the equation and, secondly, for ehminating A, 
Eq. fl22|) is equivalent to the following one: 



p U 



y_iy_ 

p D 



0. (28) 



Notice that previous equation involves only / when p and Q are written explicitly in terms 
of / from their definitions. 

We have found a decoupled system of quasi-linear PDEs fl25]) . (1271) and fl25]) according 
to the following steps: First, Eq. ( |28|) can be solved for /. Second, Eq. (!25|) can be solved 
for A. Third, Eq. ( 127|) can be solved for 6. Finally, a can be obtained from Eq. (12^ . 
Assuming that all the metric components A, B, C, D are continuously different iable func- 
tions, the initial value problem with respect to this set of equations has always local (both 
in space and time) solution.!^ We have reached the main objective of this section, that 
is to prove the following result. 
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Theorem: Any spherically symmetric spacetime can be locally sliced by a family of 
maximal spacelike hypersurfaces. 

Let us make a few comments about first order partial differential equations. If the 
functions appearing in the equation are continuously differentiable and we are looking for 
smooth solutions, given continuously differentiable initial data, it always exists a solution 
of the equation in a domain containing the given initial data. The domain of the solution 
can be the domain of initial parameters or a smaller one; so that is why, in principle, one 
can only guarantee local existence of solutions. Moreover, in this kind of equations one 
can use the so-called "characteristic method" in order to obtain, analytical or numerically, 
solutions of the equation. This method is based on two facts: i) given initial data at a 
point, it always exists a characteristic curve containing the given initial data; ii) a solution 
of the equation can be constructed by the union of a set of characteristic curves. So, a 
solution of the equation can be obtained by solving a set of characteristic curves around 
a given point. 

Let us also point out that in the resolution of the characteristic curves some func- 
tional relations involving the dependent and independent variables have to be constant; 
the cardinal of the set of independent functional relations depends on the number of in- 
dependent variables of the equation, and the general solution of the first order partial 
differential equation can be written as an implicit function of those ones. 

Once the local existence of maximal slicings has been proved, we are going (see next 
Section) to design a practical procedure to build up these slicings. 

5 A procedure to construct maximal slicings 

The proof of the above theorem was based on the explicit finding of the minimum set of 
equations to be solved in order to obtain maximal slicings in SSSTs. Let us analyze a 
constructive algorithm that complements the theorem. We consider two different cases, 
/ = and / ^ 0. 

5.1 First case / = 

If / = 0, Eq. fl28p is reduced to = 0, the maximal slicing condition provided by 
Eq. (H]). In this case, from Eqs. ([27]), (j^S]) and ({21]), it turns out that h{t) ^ arbitrary. 
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X = Z{t) a/D exp I — e J \/ B / Ddr j , with Z(t) > arbitrary and R constant, and 

R 



a = bVD 



r 

Z/Z + eC/VBD -e J {^^J B / Dy dr 



R 



The fields X and Y can be integrated according to: 

Proposition 1: If the hypersurfaces t = constant in the spacetime defined by the me- 
tric (QP are maximal, the following change of coordinates 



i = i{t), f = Z-^ exp j e J ^/bJ 



Ddr 



(29) 



R 



with t = b ^ , Z{t) > and R constant, allows one to write the metric in isotropic form 

2 



r \ B r 



eb D 



dp + 2 —\ — C - eVSD -\dtdf 
r \ B \ r 



D 



+^(^dP + Pdn^). 



(30) 



This metric form is adapted to the same initial maximal slices, that are now labeled as 
if: = constant. It remains invariant under the inversion f — )■ T(t)/r when, in addition, the 
sign of the radial component of the field Y is reversed by the change e — )■ — e. 



5.2 Second case / 7^ 

When / 7^ 0, the explicit expressions of P/f and Q in terms of /, allow us to define the 
variable 



^ f = eF 



B 



(31) 



^WTb V 1 - 

Then, using the definition of P and Q given by Eqs. ffTTl) and f|T2|) respectively, we have 



P _1-FC 
1~ FB ' 



Q 



-C + F/2 



p = e 



1 - PF^ 



B 



(32) 
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Eqs. ( 127|) . (125!) and (!28|) are rewritten, respectively, as 

1 



FB 



B 



1 - 



,2 ^ ^^^(2 -FC)-C F' 



B 



+F\l 



2\l 



B 



A + 



FB 



\' = \ 



D 1-FC D' eVl - 
2D ^ ~Tb~2D ~ fVBD 



(33) 



(34) 



2V 



FB{P)-+ {2-PF^ - FC){1 

Wb 

B BD {C - Fl'^){DB' -2BD') 



2BD 



C 



(35) 



where L = (1 - PF'^)/{fFB). The solution can be obtained numerically solving the 
quasilinear Eqs. fl35l) . flMj) and fl33l) for, respectively, F, A and 6, given initial data for 
t = to. These equations are of hyperbolic type, analogous to the well-known advective 
equation in classical fluid dynamics. The three equations have a common "advective 
velocity" : = (1 — FC) / (FB). This velocity can be interpreted as the radial component 



of the light velocity in the modified metric ds = ds -dt . 

J 

On the other hand, the coordinate 1-forms dt and dr are easily written from Eq. (ITSll . 
taking into account Eqs. fl32|) . Then, we arrive to this result: 

Proposition 2: For any given SSST a maximal slicing can be built up according to the 
following steps: 

(i) To start with the metric form (QP, with f = C"^ - AB > 0, 5 > and D > 0, to 
solve Eq. (3^) for F , and then, to solve Eq. ^3B\j for h. 



(36) 



(ii) To define 

r = 6-1(1 - fF^)-^ [(1 - FC)dt - FBdr] 
which is necessarily a closed 1-form. 

(Hi) To find a potential t{t,r) ofT, T = dt. Then, the level surfaces oft are spacelike 
and maximal. 



Analogously, by applying Eq. ( 1211) . r can be integrated once F and A are solved. The 
change of coordinates {t = i{t,r),f = f{t,r),9,ip^ is completely defined from the fields 
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X, Y given by Eqs. (fT3|) and (JH]). Then, the metric form can be deduced in the new set 
of coordinates. The new shces are maximal and Proposition 1 can again be apphed. 



6 Maximal slicings in flat spacetime 



In order to illustrate how the procedure designed in previous section works, let us consider 
the simple case of the Minkowski spacetime starting with inertial spherical coordinates, 
ds'^ = —dt"^ + dr'^ + PdVt^. According to Proposition 1, t = with t = b~^, f = r'' 
and the resulting metric is 

2" 



ds^ 



^^dZ\ 



dt J 



dp + 2 {r Z f-' ^didr + (P-I Z^) ' {dr^ + r^dQ^) , (37) 



b = b(t) ^ 0, Z = Z{t) > being arbitrary functions. So, we have the same slicing that 
the initial inertial one (this foliation is totally geodesic, i.e., its extrinsic curvature tensor 
vanishes). However, other maximal slicings (non-inertial, in general) can be obtained as 
a consequence of Proposition 2. According to it, Eq. f p5|) reads 



F + F-^F' = -2r-Hl-F'^). 



(38) 



Following the last comments of Sec. HI the general solution of this equation can be given 
in terms of an implicit function. 



$1 (m{F) + t (F-2 - l) , r'' (F'^ - l) 



1/4 



0, 



with < < 1, m{F) = 'i^/\F\ 2F1 (1/4, 3/4; 5/4; F^), 2F1 being the classical standard 
hypergeometric series^ and $1 being an arbitrary function. Once F is known, the general 
solution of Eqs. (133|) and (IMIl are, respectively, $2(^5-^2) and $3(11,13), where $2 and $3 
are arbitrary functions, and 



/i = r - y F-^dt, 


I2 = exp 



:i + F2)F 



+ 2F' 



dt 



1 -F2 



/g = A"^ exp j r-^ j F-\l - eVl - F'^)dt 



(39) 
(40) 
(41) 
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This study can be extended to the Schwarzschild spacetime using different metric forms, 

specially those that are regular at the horizon, like the Painleve-Gullstrand form and other 

similar ones.^ Maximal slincings can also be constructed in other spacetimes (Reissner- 

Nordtrom, Vaidya, Lemaitre-Tolman-Bondi, etc.) and the results compared with previous 
„.Q^.1,3 I2UI21I22I23I21 



7 Maximal slicings in closed Friedmann spacetime 

The pioneering analysis by Oppenheimer and SnydeiJS^lon the dynamics of the collapse of a 
homogeneous sphere of pressure- zero matter allowed them to obtain an analytical solution 
of Einstein equations for a dynamical spacetime. Its metric, a special case of Lemaitre- 
Tolman-Bondi metric, has a simple form (see below), being the Friedmann metric in the 
interior. We are going to apply the procedure designed in Sec. [5] to build maximal slicings 
in a closed Friedmann spacetime. We start with the following metric form: 

ds'^ = [-drf + dx^ + sin^ x^fi^) , (42) 

where R = Rq{1 — cost]), Rq being a positive constant, and r] E (— 7r,0), x ^ (0,7r), 
being the domain of the parameters {t] and x play the role of t and r in Eq. ([T]) 
respectively). The hypersurfaces rj = constant are not maximal. We are going to solve 
numerically Eqs. ([28]), ([27D and ([25]). Let us define / := Rf. 
Eq. ^ reads 

f df e df ri 2ef cot X n 

+ , ~ + Qcot^+ ; ^ = 0. (43) 



Without loss of generality, it can be considered e = 1, since for e = — 1 the solution will 
be — /. Previous equation is then equivalent to 

/ 9f ^df^ -6^/1 + /2 cot ^ - 2/ cot X- (44) 



l^p9v dx V 2 

We use the characteristic method for solving the above equation (l44l) . Let {ri{fi), x(/i), /(/i)) 

be a characteristic curve depending on a parameter /i such that — ^ = — , , = 1, 

afi w ]^ _|_ j2 ""1^ 



df I ~ 7] 

and — = — 6\/ 1 + cot 2/ cot x- Substituting dx = da, and considering then 

djji ^ 2 
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Figure 1: Projection onto the t] ~ x plane of characteristic curves of Eqs. (H5il . with 
fix = 7r/2) = 0, ■r]{x = 7r/2) = -kn, < k < 1, initial data. 



(?7(x), fix)): the characteristic system associated to Eq. (144|) is written as 



dv/dx = //Vl + df/dx = -6Vl + P cot | - 2/ cot 



(45) 



Notice that, although rj = —n is out of the range of the parameter, f = 0, r] = — vr is a 
solution of the previous system. Notice also that solutions of Eqs. ( 115]) satisfy f (^'^ ~ ^ 
—f ^— + and T] — = r] + x^. Eqs. f HSj) can be solved numerically, imposing 



/(X = 7r/2) = 0, r]ix = v^/2) = -A; vr, < A; < 1, 



(46) 



as initial data. In Fig. [T]we plot the projection onto the 1] — x plane of the characteristic 
curves, for the values k = —1/50, 1 = 1,..., 50. In Fig. [2] we plot the projection onto the 
f — X plane of the characteristic curves, for the values k = 1/5, I = 1, ... ,4. The union 
of the characteristic curves, for all the possible values of fc, < A; < 1, covers the whole 
domain of parameters for t] and x- This union constitutes a global solution of Eq. ( H4|) 
for /. 

Once / (and therefore /) is solved, we focus on Eq. ( |271) for b, which reads (after 
multiplying by eR) 



ef 



~ (9 log 6 



df] 



^ ^ p d^ogh ^ df_ 



f 



+ e 1 



+ J ^. (47) 



/ 



+ P 

From Eq. (1441) . the combinations ef and will remain the same for both values of 
e. Without loss of generality, we take e = 1 and the previous solution for /. Due to 
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-2.xl0'h ' 
-4.X10' 
-6.X10- 



I 4 



n5 . 



3ir 
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Figure 2: Projection onto the f — x plane of characteristic curves of Eqs. ( H5ll . with 
f(x = n/2) = 0, r/(x = n/2) = -kn, k = 1/5 (sohd hue), 2/5 (dashed hue), 3/5 
(dash-dotted hue), 4/5 (dotted hne) initial data. 

df . df . 

numerical reasons, we substitute (using Eq. (l44l) ) — in terms of — , leading to 

ox or] 

f d log b ^ d log b 



1 + / 



2 dr] 



dx 




f- + P 

We use again the characteristic method for solving the above equation ( H8 
consider a characteristic curve (?7(x), log6(x)), such that 

dr]/dx 



(48) 
Let us 



dlogb 
dx 




f 



+ 1 



2f 



rT cot X - 



6/ 



cot-. 



(49) 



Notice that the projection onto the — x plane will be exactly the same as the one for the 
characteristic curves associated to /. In practice, the characteristic curves for log 6 can 
be calculated simultaneously to the ones for /. Notice also that, although r] = — vr is out 
of the range of the parameter, b = 1, t] = — vr (with / = 0) is a solution of the previous 
system. Eqs. (HH]) can be solved numerically, imposing 

log6(x = vr/2) = 0, (50) 

and the previous values for / and rj as initial data. In Fig. |3]we plot the projection onto 
the log 6 — X plane of the characteristic curves, for the values k = 1/5, I = 1, ... ,4. The 
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Figure 3: Projection onto the log 6 — x plane of characteristic curves of Eqs. (H9ll . with 
log6(x = 7r/2) = 0, previous values for / and r], k = 1/5 (solid line), 2/5 (dashed line), 
3/5 (dash-dotted line), 4/5 (dotted line), initial data. 



union of the characteristic curves, for all the possible values of k,0 < k < 1, covers the 
whole domain of parameters for t] and x- This union constitutes a global solution of 
Eq. (iH]) for log b. b ^ is satisfied. 



We focus now on ( |25|) for A, which reads (after multiplying by eR/ y 1 + P) 



ef dhgX Slog A 2ef r] cosx-t/dl + p 

H = — , cot - H : — 

2 sm X 



The combinations ef and will remain the same for both values of e, so we can use 
previous solution for / but we have to keep e in the last term. 



(51) 



/ dlogX Slog A 2/ 7] cosx-e/Vl + /^ 
H 7^ = — , cot - H r 



drj 



dx 



(52) 



+ -'I -A ^i + /2 ^ smx 

Therefore, we will have two solutions log A^ which do not differ only in the sign. Using 
again the characteristic method for solving the above equation fl52|) . we consider a cha- 
racteristic curve (//(x), log A(x)), such that 



dv/dx = f/\^l + f 



^2 dlogX _ 2/ 

^ p 2 ' sm X 



dx 



ri cosx - e/-y/l + /2 
cot- + - 



(53) 



The projection onto the f] — x plane is exactly the same as the one for the other charac- 
teristic curves. Hence, the corresponding projections for log A can be calculated simul- 
taneously to them. Notice that, although rj = —n is out of the range of the parameter. 
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(a) e = 1. (b) e = -1. 

Figure 4: Projection onto the log A — x plane of characteristic curves of Eqs. f l53|) . with 
logA(x = 7r/2) = 0, previous values for / and r/, k=l/b (solid line), 2/5 (dashed line), 
3/5 (dash-dotted line), 4/5 (dotted line), initial data. 

Ae = Ao sin X (cot x/2)% rj = — vr, Aq being a positive constant, is a solution of the previous 
system (with / = 0). Eqs. fl53|) can be solved numerically, imposing 



logA(x = 7r/2)=0, 



(54) 



and the previous values for / and rj as initial data. These initial data corresponds to the 
value Aq = 1 in the limit rj = —it. In Fig. H] we plot the projection onto the log A — x plane 
of the characteristic curves, for the values k = 1/5, I = 1, ... ,4, and for both values of e. 
The union of the characteristic curves, for all the possible values of k,0 < k < 1, covers 
the whole domain of parameters for t] and x- This union constitutes a global solution of 
Eq. ([52]) for log A. A > is satisfied. 

Notice that the fact that the projection onto the 1] — x plane of the characteristic 
curves for /, log 6 and log A coincides is just a consequence of the fact that Eqs. (I551) - fl5^ 
share the same advective velocity. 

Once /, b and A have been derived, we can calculate the unknown a from Eq. (!24|) . 

bR 



sinx {dlogX/drj — 2cotri/2) — f 



1 + P 

Also, from Eq. (^T^, we can directly obtain 

. R . 

r = — sm Y , 
A 



(55) 



(56) 



In Fig. |5]we plot the projection onto the r — x plane of the characteristic curves, for the 
values k = 1/5, I = 1, ... ,4, and for both values of e. 
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X 



X 



(c) A; = 3/5. (d) = 4/5. 

Figure 5: Projection of characteristic curves onto f — x plane. Solid lines correspond to 
e — 1, while dashed hues correspond to e = — 1. 
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10 




Figure 6: Projection onto the t — x plane of characteristic curves, with t{x = 7r/2) = 
1], rj{x = ^/2) = —kn, k = 1/5 (sohd hne), 2/5 (dashed hne), 3/5 (dash-dotted hne), 4/5 
(dotted hne) initial data. 



The conditions for the variable t are the ones in Proposition 2, once F is written in 
terms of / and the components of the initial metric fH2l) in this particular case. From 
Eq. ( 136|) (let us draw reader's attention to the fact that we are using t] and x instead of t 
and r, respectively), and the corresponding values of / and b, i has to satisfy di/dx\ 



I rj=—iT 



0, i.e., t\ 



has to be constant. Moreover, dt/dr]\ 



1. We can choose consistently 



T] and t\ = —71. The values of t in the whole domain of parameters can 
be numerically obtained from the expression of di/drj, coefficient of dt] in Eq. (1361) for F, 
and the condition 1 1 _ ,„ = r]- In Fig. [6] we plot the projection onto the i — x plane of 
the characteristic curves, for the values t{x = 7r/2) = t], r]{x = tt/2) = —kn, k = 1/5, 

1 = 1,. ..A- 

Numerically, the metric ( H2ll can be written in the new set of coordinates ^t,r,9,ip^, 
evaluating the corresponding Jacobian matrix from (fT5|) . The hypersurfaces t = constant 
are now maximal. 



8 Conclusions 

We have proved that any SSST can be locally sliced by a family of maximal spacelike 
hypersurfaces. 

We have designed a general geometrical method to build up maximal slices in any 
SSST by solving three decoupled first order quasi-linear PDEs fl25l) . fl27|) and fl28|) . It has 
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been applied in the cases of Minkowski (vacuum) and Friedmann (matter) spacetimes. In 
the first case, a flat spacetime, the equations leading to obtain maximal slicings, which are 
not inertial in general, can be solved analytically. In the second case, closed Friedmann 
spacetime, the corresponding equations leading to obtain maximal slicings are solved 
numerically. 

An interesting by-product for Numerical Relativity of the approach presented in this 
paper has to do with the assessment of 3D codes written, as customary, in Cartesian 
coordinates. Let us consider two codes NCI and NC2 such that only NCI uses a gauge 
which is maximal. Hence, the evolution with code NC2 of any initial data admitting a 
spherically symmetric limit could be compared to the evolution produced by code NCI, by 
simply using our procedure to generate a SSST satisfying the maximal slicing condition. 
To emphasize this point we have applied our procedure to build up, numerically, maximal 
slicings in a Friedmann spacetime. If one takes into account that Oppenheimer-Snyder 
collapse is one of the basic test-beds for any numerical code aiming to solve Einstein 
equations, our analysis paves the way to calibrate numerical codes that use formalisms 
that do not incorporate the maximal slicing condition. 

We have outlined a method that should be extended to deal with the construction of 
constant mean curvature slicings. The results on maximal slicings we have presented here 
can be considered as a step to gain some insights in this direction. 

Let us make a comment about maximal hypersurfaces in a cosmological scenario. In 
the seventies, Brill,!^ and Brill and FlahertjJ^ analyzed the existence and uniqueness 
of maximal hypersurfaces in a cosmological context. Brill concludes that closed surface 
of stationary area (MS) in closed universes are exceptional and scarce^^ in particular. 
Brill finds: Any closed MS in a closed universe which satisfies Einstein's equations with the 
energy condition is a local maximal surface and is locally isolated^^These authors reached 
these conclusions in analyzing non-static Friedmannian cosmologies, for which such a 
hypersurfaces are identified with instants of cosmic time. In fact, in these cosmologies, 
the volume functional (relative to the matter flow) might take a critical value only if the 
expansion of the cosmological observers vanishes, i.e. at a particular cosmological moment, 
depending on the particular closed universe. The above conclusions by Brill and Flaherty 
and the one presented in this paper, i.e., the local existence of a maximal slicing in any 
SSST, are compatible. In fact, the result by Bill and Flaherty involves the hypothesis of 
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global spatial compacteness meanwhile our result concerns only local existence, and no 
global topological property has been required. 

As we have pointed out from the very beginning, constructing maximal slicing are 
closely related to find a vorticity-free observer without expansion. Such an Eulerian ob- 
server is not necessary to be geodesic or shear- free, i.e., it is not (generically) a cosmological 
observer moving with the matter content of the universe. Then, as a corollary of the re- 
sults we have presented, maximal slicings exist locally (and can be constructed) in any 
Robertson- Walker spacetime: they are those shcings which are associated with Eulerian 
observers without expansion. 

Our statement about the local existence of a maximal slicing in any SSST, and the 
examples of constructions we have considered, support the physical interest in considering 
such a local foliations when their global existence is forbidden from topological require- 
ments. 

Acknowledgements 

I. C.-C. acknowledges support from the Alexander von Humboldt Foundation. We wish to 
thank Dr. Jose Luis Jaramillo for detailed appreciations and interesting comments about 
this paper and related issues. This work has been also supported by the MICINN Grant 
No. AYA2007-67626-C03-01 and the MICINN-FEDER project No. FIS2009-07705. 

References 

^ A. Lichnerowicz, "L'integration des equations de la gravitation relativiste et le 
probleme des N corps" J. Math. Pure Appl. 23 37-63 (1944). 

^ Y. Choquet-Bruhat, "Maximal submanifolds and sub manifolds with constant mean 
extrinsic curvature of a lorentzian manifold," Annali della Scuola Normale Superi- 
ore di Pisa Classe di Scienze 4^ serie, tome 3, number 3, 361-376 (1976). 

^ M. Cantor, A. Fisher, J. Marsden, N. O Murchadha, and J. M. York, "The existence 
of Maximal Slicings in Asyptotically Flat Spacetimes," Commun. Math. Phys. 49, 
187-190 (1976). 

^ J. E. Marsden, and F. J. Tipler "Maximal hypersurfaces and foliations of constant 



22 



mean curvature in general relativity," Pliys. Rep. 66, 109-139 (1980). 

^ R. Bartnik, "Existence of maximal surfaces in asymptotically flat spacetimes," 
Commun. Math. Phys. 94, 155-175 (1984). 

^ P. T. Chrusciel, and R. W. Wald, "Maximal hypersurfaces in stationary asymptot- 
ically fiat spacetimes," Commun. Math. Phys. 163, 561-604 (1994). 
G. A. Burnett, and A. L. Rendall, "Existence of maximal hypersurfaces in some 
spherically symmetric spacetimes," Classical and Quantum Gravity, 13, 111-123 
(1996). 

^ Y. Choquet-Bruhat, General Relativity and the Einstein Equations (Oxford Uni- 
versity Press, Oxford, 2008). 

^ L. Smarr, and J. W. York, "Kinematical conditions in the construction of space- 
time," Phys. Rev. D 17, 2529-2551 (1978). 

L. Smarr, and J. W. York, "Radiation gauge in general relativity," Phys. Rev. D 
17, 1945-1956 (1978). 

M. Shibata, and T. Nakamura, "Evolution of three-dimensional gravitational 
waves: Harmonic slicing case," Phys. Rev. D 52, 5428-5444 (1995). 
J. A. Isenberg, "Waveless Approximation Theories of Gravity," Int. J. Mod. Phys. 
D 17, 265-273 (2008). 

S. Bonazzola, E. Gourgoulhon, P. Grandclement, and J. Novak, "Constrained 
scheme for the Einstein equations based on the Dirac gauge and spherical co- 
ordinates," Phys. Rev. D, 70 104007 (2004). 

L-M. Lin, and J. Novak, "Rotating star initial data for a constrained scheme in 

numerical relativity," Classical and Quantum Gravity 23, 4545-4561 (2006). 

E. Gourgoulhon, 3+1 Formalism and bases of Numerical Relativity, Lectures at 



Institute Henri Poincare (Paris, 2006), arXiv:gr-qc/0703035 

M. Alcubierre, Introduction to 3+1 Numerical Relativity (Oxford University Press, 
Oxford, 2008). 

R. Schoen, and S. T. Yau, "On the proof of the positive mass conjecture in general 
relativity," Commun. Math. Phys. 65, 45-76 (1979). 

R. Schoen, and S. T. Yau, "Proof of the positive mass theorem. II." Commun. 
Math. Phys. 79, 231-260 (1981). 

J. W. York, "Kinematics and dynamics of General relativity," Sources of Gravita- 



23 



tional Radiation, Edited by L. L. Smarr, (Cambridge University Press, Cambridge 
1979) pp. 83-126. 

F. Estabrook, H. Wahlquist, S. Christensen, B. DeWitt, L. Smarr, and E. Tsiang, 
"Maximally Slicing a Black Hole," Phys. Rev. D 7, 2814-2817 (1973). 
R. Beig, and N. O Murchadha, "Late time behavior of the maximal slicing of the 
Schwarzschild black hole," Phys. Rev. D 57, 4728-4737 (1998). 

B. L. Reinhart, "Maximal foliations of extended Schwarzschild space," J. Math. 
Phys. 14, 719 (1973). 

E. Malec, "Event horizons and apparent horizons in spherically symmetric geome- 
tries," Phys. Rev. D 49, 6475-6483 (1994). 

L. I. Petrich, S. L. Shapiro, and S. A. Teukolsky, "Oppenheimer-Snyder collapse 
with maximal time slicing and isotropic coordinates," Phys. Rev. D 31, 2459-2469 
(1985). 

D. M. Eardley, and L. Smarr, "Time functions in numerical relativity: Marginally 
bound dust collapse," Phys. Rev. D 19, 2239-2259 (1979). 

J. Guven, and N. O Murchadha, "Flat fohations of spherically symmetric geome- 
tries," Phys. Rev. D 60, 104015 (1999). 

V. Husain, A. Qadir, and A. A. Siddiqui, "Note on flat foliations of spherically 

symmetric spacetimes," Phys. Rev. D 65, 027501 (2001). 

R. Beig, and A. A. Siddiqui, "Uniqueness of flat spherically symmetric spacelike 
hypersurfaces admitted by spherically symmetric static spacetimes," Classical and 
Quantum Gravity 24, 5435-5439 (2007). 

C. Doran, "New form of the Kerr solution," Phys. Rev. D 61, 067503 (2000). 

J. Natario, "Painleve-Gullstrand coordinates for the Kerr solution," Gen. Relativ. 
Gravit. 41, 2579-2586 (2009). 

C-Y. Lin and C Soo, "Generahzed Painleve-Gullstrand metrics," Phys. Lett. B 
671, 493-495 (2009). 

A. Herrero, and J. A. Morales-Lladosa, "Painleve-Gullstrand synchronizations in 
spherical symmetry," Classical and Quantum Gravity 27, 175007 (2010). 

E. Malec, and N. O Murchadha, "General spherically symmetric constant mean 
curvature foliations of the Schwarzschild solution," Phys. Rev. D 80, 024017 
(2009). 



24 



J. R. Oppenheimer, and H. Snyder, "On Continued Gravitational Contraction," 
Phys. Rev. 56, 455-459 (1939). 

J. J. Ferrando, and J. A. Saez, "An intrinsic characterization of spherically sym- 
metric spacetimes," Classical and Quantum Gravity 27, 205024 (2010). 
J. Plebanski, and A. Krasinski, An Introduction to General Relativity and Cosmol- 
ogy (Cambridge University Press, Cambridge, 2006). 

L. I. Petrich, S. L. Shapiro, and S. A., Teukolsky, "Oppenheimer-Snyder collapse 
in polar time slicing," Phys. Rev. D 33, 2100-2110 (1986). 

E. Gourgoulhon, Ph. D. Thesis, Universite de Paris Vll (1992). 

J. V. Romero, J. M. Ibaficz, J. M. Marti, and J. A. Miralles, "A New Spherically 

Symmetric General Relativistic Hydro dynamical Code," Astrophys. J. 462, 839 

(1996). 

I. Cordero-Carrion, J. M. Ibanez, and J. A. Morales-Lladosa, "On the local exis- 
tence of maximal slicings in spherically symmetric spacetimes," J. Phys.: Conf. 
Ser. 228, 012055 (2010). 

E. Zauderer, Partial differential equations of applied mathematics, Second edition. 
(John Wiley & Sons, New York, 1989). 

M. Abramowitz, and I. A. Stegun, Handbook of mathematical functions, p. 556, 
(Dover Publications, Inc., New York, 1965). 

P. Papadopoulos, and J. A. Font, "Relativistic hydrodynamics around black holes 
and horizon adapted coordinate systems," Phys. Rev. D 58, 024005 (1998). 
D. Brill, "Maximal surfaces in closed and open spacetimes," in Proceedings of the 
1st Marcel Grossmann meeting on general relativity. Edited by R. Ruffini (North- 
Holland Pub. Comp., Amsterdam, 1977) pp. 193-206. 

D. Brill, and F. Flaherty, "Isolated Maximal surfaces in Spacetime," Commun. 
Math. Phys. 50, 157-165 (1976). 



25 



